Excitation Spectrum Gap and Spin- Wave Velocity of XXZ Heisenberg Chains: 
Global Renormalizat ion- Group Calculation 
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The anisotropic XXZ spin-^ Heisenberg chain is studied using renormalization-group theory. The 
specific heats and nearest- neighbor spin-spin correlations are calculated thoughout the entire tem- 
perature and anisotropy ranges in both ferromagnetic and antiferromagnetic regions, obtaining a 
global description and quantitative results. We obtain, for all anisotropics, the antiferromagnetic 
spin-liquid spin-wave velocity and the Isinglike ferromagnetic excitation spectrum gap, exhibiting 
the spin-wave to spinon crossover. A number of characteristics of purely quantum nature are found: 
The in-plane interaction sfsj + s v s y induces an antiferromagnetic correlation in the out-of-plane sf 
component, at higher temperatures in the antiferromagnetic XXZ chain, dominantly at low temper- 
atures in the ferromagnetic XXZ chain, and, in-between, at all temperatures in the XY chain. We 
find that the converse effect also occurs in the antiferromagnetic XXZ chain: an antiferromagnetic 
sfsj interaction induces a correlation in the sf H component. As another purely quantum effect, 
(i) in the antiferromagnet, the value of the specific heat peak is insensitive to anisotropy and the 
temperature of the specific heat peak decreases from the isotropic (Heisenberg) with introduction 
of either type (Ising or XY) anisotropy; (ii) in complete contrast, in the ferromagnet, the value and 
temperature of the specific heat peak increase with either type of anisotropy. 

PACS numbers: 67.40.Db, 75.10.Pq, 64.60.Cn, 05.10.Cc 



I. INTRODUCTION 

The quantum Heisenberg chain, including the possi- 
bility of spin-space anisotropy, is the simplest nontrivial 
quantum spin system and has thus been widely studied 
since the very beginning of the spin concept in quan- 
tum mechanic s El 13 , H| . Interest in this model continued 
d, H, H, 0, H, l9l. Ildf and redoubled with the exposition 
of its richly varied low- temperature behavior [ill. [l2l [l3j 
and of its relevance to high-temperature superconductiv- 
ity El, [H, E3, Ell • It has become clear that antifer- 
romagnetism and superconductivity are firmly related to 
each other, adjoining and overlapping each other. 

A large variety of theoretical tools have been employed 
in the study of the various isotropic and anisotropic 
regimes of the quantum Heisenberg chain, including 
finite-systems extrapolation @, El] , linked-cluster 0] and 
dimer-cluster [2(| expansions, quantum decimation [2l| . 
decoupled Green's functions 22], quantum transfer ma- 
trix [23, [24[ , high-temperature series expansion [25| , and 
numerical evaluation of multiple integrals [26| . The XXZ 
Heisenberg chain retains high current interest as a theo- 
retical model [13, HH with direct experimental relevance 

In the present paper, a position-space renormalization- 
group method introduced by Suzuki and Takano [13, HH 
for d = 2 dimensions and already applied to a number 
of d > 1 systems [13, H [H \M, H HE HI is used 
to compute the spin-spin correlations and the specific 
heat of the d = 1 anisotropic quantum XXZ Heisenberg 
model, easily resulting in a global description and de- 
tailed quantitative information for the entire temperature 
and anisotropy ranges including the ferromagnetic and 



antiferromagnetic, the spin-liquid and Isinglike regions. 
We obtain, for all anisotropics, the antiferromagnetic 
spin- liquid spin- wave velocity and the Isinglike ferromag- 
netic excitation spectrum gap, exhibiting the spin- wave 
to spinon crossover. A number of other characteristics of 
purely quantum nature are found: The in-plane interac- 
tion sfs^ + sfs^ induces an antiferromagnetic correlation 
in the out-of-plane sf component, at higher temperatures 
in the antiferromagnetic XXZ chain, dominantly at low 
temperatures in the ferromagnetic XXZ chain, and, in- 
between, at all temperatures in the XY chain. We find 
that the converse effect also occurs in the antiferromag- 
netic XXZ chain: an antiferromagnetic sfs| interaction 
induces a correlation in the sf" component. As another 
purely quantum effect, (i) in the antiferromagnet, the 
value of the specific heat peak is insensitive to anisotropy 
and the temperature of the specific heat peak decreases 
from the isotropic (Heisenberg) with introduction of ei- 
ther type (Ising or XY) anisotropy; (ii) in complete con- 
trast, in the ferromagnet, the value of the specific heat 
peak is strongly dependent on anisotropy and the tem- 
perature of the specific heat peak increases with either 
type of anisotropy. This purely quantum effect is a pre- 
cursor to different phase transition temperatures in three 
dimensions [H, US H3, [HI- Our calculational method 
is relatively simple, readily yields global results, and is 
overall quantitatively successful. 
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II. THE ANISOTROPIC QUANTUM 
HEISENBERG MODEL AND THE 
RENORMALIZATION-GROUP METHOD 

A. The Anisotropic Quantum Heisenberg Model 

The spin-i anisotropic Heisenberg model (XXZ model) 
is defined by the dimensionlcss Hamiltoiiian 
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TABLE I: The two-site basis eigenstates that appear in 
Eq.©. These are the well-known singlet and triplet states. 
The state \(j>3) is obtained by spin reversal from \<f)i), with the 
same eigenvalue. 



-Pi = { [ J *y + + J * s i s T\ + G ) . (!) 

where (3 = 1/ksT and (ij) denotes summation over 
nearest-neighbor pairs of sites. Here the s™ are the quan- 
tum mechanical Pauli spin operators at site i. The ad- 
ditive constant G is generated by the renormalization- 
group transformation and is used in the calculation of 
thermodynamic functions. The anisotropy coefficient is 
R = J z /J xy . The model reduces to the isotropic Heisen- 
berg model (XXX model) for \R\ = 1, to the XY model 
for R = 0, and to the Ising model for \R\ — > oo. 



e -p'H'(i,k) = rp,. e {-/m{i,j)-fm(j,k)} } (4) 

where i,j, k are three successive sites. The operator 
— /3'7i'(i,fc) acts on two-site states, while the operator 
—(3H(i,j) — flH-ij, k) acts on three-site states, so that we 
can rewrite Eq.(j4]) in the matrix form, 



[UiV k e 



Ui wj v k \e 



-f3H(t,j)-pn(j,k)\ 



H w 3 v kl > 



(5) 



B. Renormalization-Group Recursion Relations 

The Hamiltonian in Eq.((T]) can be rewritten as 



■pH = '£{-pH(i,i + i)}. 



(2) 



where j3TL{i, i + 1) is a Hamiltonian involving sites i and 
i + 1 only. The renormalization-group procedure, which 
eliminates half of the degrees of freedom and keeps the 
partition function unchanged, is done approximately [13, 



Tr odd e"^ =Tr odd e^-^' l+1 )> (3) 
=Tr e Y.° Ad {-Pn(i-i,i)-t3H{i,i+i)} 



, ^ r rr , e {-f3H(,i-U)-0K{i,i+l)} 



A e -/3'w'(i-i,i+i) _ e E° d,1 {-/5'H'(.-i,.+i)} = e -P'n\ 



Here and throughout this paper, the primes are used for 
the renormalized system. Thus, at each successive length 
scale, we ignore the non-commutativity of the operators 
beyond three consecutive sites, in the two steps indicated 
by ~ in the above equation. Since the approximations 
are applied in opposite directions, one can expect some 
mutual compensation. Earlier studies [13, H2, HE H3j 
Hf| have been successful in obtaining finite-temperature 
behavior on a variety of quantum systems. 
The transformation above is summarized by 



where state variables u, v, w, u, v can take spin-up or spin- 
down values at each site. The unrenormalized 8x8 ma- 
trix on the right-hand side is contracted into the renor- 
malized 4x4 matrix on the left-hand side of Eq. ([5| . We 
use two-site basis states vectors {|</> P )} and three-site ba- 
sis states vectors {IV'g)} to diagonalize the matrices in 
Eq.®. The states {|</> p )}, given in Table HI are eigen- 
states of parity, total spin magnitude, and total spin z- 
component. These {|0 P )} diagonalize the renormalized 
matrix, with eigenvalues 



Ai — -Jz + G', A 2 — +-Jxy - jJz + G' , 



(6) 



The states {|V>g)}j given in Table HT1 are eigenstates of 
parity and total spin z-component. The {IV'p)} diagonal- 
ize the unrenormalized matrix, with eigenvalues 



Ai — -J z + 2G, A4 — 2G, 



(7) 

A 2 = -- (j z + v /8^ + J?) + 2G, 
With these eigenstates, Eq.© is rewritten as 

lp = (4>p\e~ ' U ' (i,k) \(l>p) = ^p\ U i V k){UiWjVk\^q)- 

v,w,q 

(^le-^WJ-^W'*^,)^,^^)^*!^,) . (8) 
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TABLE II: The three-site basis eigenstates that appear in 
Eq.JHJ) with coefficients a = {-J z + ^/8J| H + J'i)/2J xy ,T = 
(J z + ^d>J% v + Ji)/2J xy and normalization factors /j,, v. The 
states | i/jg_s } are obtained by spin reversal from |t/>i_4), with 
the same respective eigenvalues. 



renormalization-group trajectories do not change. In 
fact, this is part of a more general symmetry of the XYZ 
model, where flipping the signs of any two interactions 
leaves the spectrum unchanged 0]. Therefore, with no 
loss of generality, we take J xy > 0. Independent of the 
sign of J xy , J z > gives the ferromagnetic model and 
J z < gives the antiferromagnetic model. 



C. Calculation of Densities and Response 
Functions by the Recursion-Matrix Method 



Thus, there are three independent 7 P that determine the 
renormalized Hamiltonian and, therefore, three renor- 
malized interactions in the Hamiltonian closed under 
renormalization-group transformation, Eq.([T]). These 7 P 
are 



Just as the interaction constants of two consecutive 
points along the renormalization-group trajectory are re- 
lated by the recursion relations, the densities are con- 
nected by a recursion matrix T, which is composed of 
derivatives of the recursion relations. For our Hamilto- 
nian, the recursion matrix and density vector M are 



7i = e^ J * +G ' = e 2G_ H 



COsh -J8J2 +J2 



xy 



& Jxy + Jz 



j 2 = e h J £y-i J £+ G =2e 2G ~i J * 



cosh( ^8J* V + J z 2 



J z sinh 



+ 



&Jxy + Jz 



7 4 = e ~3-4VH'+ G ' = 2e 2G 



(9) 



which yield the recursion relations 
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,G' 



In (717274) 



(10) 



As expected, J xy and J' z are independent of the additive 
constant G and the derivative dcG' — b d = 2, where b = 
2 is the rescaling factor and d = 1 is the dimensionality 
of the lattice. 

For J xy = J z , the recursion relations reduce to the spin- 
| isotropic Heisenberg (XXX) model recursion relations, 
while for J xy — they reduce to spin-i Ising model re- 
cursion relations. The J z = subspace (XY model) is 
not (and need not be) closed under these recursion rela- 
tions [30I |3~H : The renormalization-group transformation 
induces a positive J z value, but the spin-space easy-plane 
aspect is maintained. 

In addition, there exists a mirror symmetry along 
the J z -axis, so that J xy {—J xy ,J z ) = J xy (J xy , J z ) and 
J z (—Jxy,J z ) — J z (Jxy,J z )- The thermodynamics of 
the system remains unchanged under flipping the in- 
teractions of the x and y spin components, since the 



T = 



dG 






dG' 

dJyy 

dJ*y 



V ® d.£ y dT I 



dG' 
dJ, 



M=(l 2( S f S f) ( S f S |>) 



(11) 



These are densities M a associated with each interaction 

K a , 



1 d\nZ 

iVa dK a 



(12) 



where N a is the number of a-type interactions and Z 
is the partition function for the system, which can be 
expressed both via the unrenormalized interaction con- 
stants as Z(K) or via the renormalized interaction con- 
stants as Z(K'). By using these two equivalent forms, 
one can formulate the density recursion relation [3^ | 







T, 



Pa 



N a dK a 



(13) 



Since the interaction constants, under renormalization- 
group transformation, stay the same at fixed points such 
as critical fixed points or sinks, the above Eq. lfTB")) takes 
the form of a solvable eigenvalue equation, 



b d M* 



M* -T. 



(14) 



at fixed points, where M — M' — M* . The fixed point 
densities are the components of the left eigenvector of the 
recursion matrix with left eigenvalue b d [391 ] . At ordinary 
points, Eq. (|13[) is iterated until a sink point is reached 
under successive renormalization-group transformations. 
In algebraic form, this means 
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(15) 



where the upper indices indicate the number of iteration 
(transformation), with M< n ) ~ M*. 

This method is applied on our model Hamiltonian. 
The sink of the system is at infinite temperature J* y = 
J* = for all initial conditions (J X y, Jz)- 

Response functions are calculated by differentiation 
of densities. For example, the internal energy is U — 
-2 (sf af ) -R^fs*), employing T = 1/J xy , and U = 
-2 (s* v s* v ) /R-(sfs*j), employing T = l/\ J z | . The spe- 
cific heat C = 8tU follows from the chain rule, 



■'•'I 



c = j: 



hi. 



3J X y 

d{2{s**sf)/R+{s^)) 



d\M 



for T = l/J xy , 

, fiorT=l/|Jt|. 

(16) 



CORRELATIONS SCANNED WITH 
RESPECT TO ANISOTROPY 



The ground-state and excitation prop erties of the XXZ 
model offer a variety of behaviors [U [H, HO, Sl| : The 
antiferromagnetic model with R < — 1 is Isinglike and 
the ground state has Neel long-range order along the z 
spin component with a gap in the excitation spectrum. 
For — 1 < R < 1, the system is a "spin liquid", with a 
gapless spectrum and power-law decay of correlations at 
zero temperature. The ferromagnetic model with R > 1 
is also Isinglike, the ground state is ferromagnetic along 
the z spin component, with an excitation gap. 

Our calculated (sfs|) and (sf* ' s^ v ) = (sf s j) = 
(sfsj) nearest-neighbor spin-spin correlations for the 
whole range of the anisotropy coefficient R are shown 
in FigtU for various temperatures. The xy correlation is 
always non-negative. Recall that we use J xy > with 
no loss of generality. In the Isinglike antiferromagnetic 
(R < —1) region, the z correlation is expectedly anti- 
ferromagnetic. As the (sfs|) correlation saturates for 
large \R\, the transverse {s^ s* v ) correlation is some- 
what depleted. In the Isinglike ferromagnetic (R > 1) 
region, the (sf s|) correlation is ferromagnetic, saturates 
quickly as the (s^ y s^ a ) correlation quickly goes to zero. 
In the spin-liquid (\R\ < 1) region, the (sf s|) correlation 
monotonically passes through zero in the feromagnetic 
side, while the (s^ y Sj V ') correlation is maximal. The re- 
markable quantum behavior of (sfs|) around R = is 
discussed in Sec.V below. It is seen in the figure that 
these changeovers are increasingly sharp as temperature 
is decreased and, at zero temperature, become discon- 
tinuous at R = 1. As seen in Fig. 1(b), at zero temper- 
ature, our calculated (sfs|) and (s^s^) correlations 
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FIG. 1: (a) Calculated nearest-neighbor spin-spin correlations 
{siSj} (thick curves from lower left) and (s^s^) (thin curves 
from upper left) as a function of anisotropy coefficient R for 
temperatures 1/J xy = 0,0.1,0.2,0.4,0.8. (b) Calcufated zero- 
temperature nearest-neighbor spin-spin correiations (thin and 
thick curves, as in the upper panel) compared with the exact 
points of Ref.0, S3, S3, EH, |3| shown with filled and open 
symbols for (sfs|) and {s* y s* v ) respectively. At R = 1, the 
calculated (sf s| \ discontinuously goes from antiferromagnetic 
to the exact result of 0.25 [40( of saturated ferromagnetism 
and the calculated (s* y s* y ~) discontinuously goes from ferro- 
magnetic to the exact result of constant zero [401 ]. 



show very good agreement with the known exact points 
0, S EE S3- Also, our results for R > 1 fully overlap 
the exact results of (sfs|) = 0.25 and « y sf ) = [If. 
We also note that zero-temperature is the limit in which 
our approximation is at its worst. 
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FIG. 2: Calculated nearest-neighbor spin-spin corre- 
lations (s^s^) (upper panels) and ( s i s j) (lower 
panels) for the antiferromagnetic XXZ chain, as a 
function of temperature, for anisotropy coefficients 
R = 0, -0.25, -0.50, -0.75, -1, -2, -4, -8, -co spanning the 
spin-liquid (left panels) and Isinglike (right panels) regions. 
Note that, in every one of the panels, the correlation curves 
cross each other. This remarkable quantum phenomenon is 
discussed in the text. 



IV. ANTIFERROMAGNETIC XXZ CHAIN 

For the antiferromagnetic XXZ chain, our calculated 
(sfs|) and (s^s^) nearest-neighbor spin-spin correla- 
tions as a function of temperature are shown in Fig[2] 
for various anisotropy coefficients R. We find that when 
J xy is the dominant interaction (spin liquid), the corre- 
lations are weakly dependent on anisotropy R. When J z 
is the dominant interaction (Isinglike), the correlations 
are weakly dependent on anisotropy R only at the higher 
temperatures. Our results are compared with multiple- 
integral results 26] in FigJ3] 

In every one of the panels of Fig[21 the correlation 
curves cross each other, revealing a remarkable quantum 
phenomenon. In a classical system, the correlation be- 
tween a given spin component (e.g., (s^sj 2 ')) is expected 
to decrease when the coupling of another spin component 
[e.g., \J Z \) is increased. It is found from the antiferromag- 
netic XXZ chain in Fig[2] that the opposite may occur 
in a quantum system: In this figure, an increase in J xy 
causes an increase in | (sfs|) | for 1/| J z \ > 0.9 and 0.4 
in the spin-liquid and Isinglike regions respectively. Con- 
versely, an increase in | J z \ causes an increase in \s^ v Sj V ) 
for l/J X y > 0.4 and 2.1 in the spin-liquid and Isinglike 
regions respectively. This quantum effect can be called 
cross-component spin correlation. 




FIG. 3: Comparison of our results (thick lines) for the cor- 
relation functions of the antiferromagnetic XXZ chain, with 
the multiple-integral results of Ref . [261] (thin lines), for var- 
ious anisotropy coefficients R spanning the spin-liquid and 
Isinglike regions. 



The antiferromagnetic specific heats calculated with 
Eq. (fTT)]) are shown in FigU] for various anisotropy coeffi- 
cients and compared, in Figs[5l [6l with finite-lattice ex- 
pansion [f| [l9| , quantum decimation plj , transfer matrix 
[24| . high-temperature series expansion [25[ results and, 
for the R — case, namely the XY model, with the exact 
result [5] C = (1/4ttT) (cos u/ cosh (^)) 2 du. The 
C(T) peak temperature is highest for the isotropic case 
(Heisenberg) and decreases with anisotropy increasing in 
either direction (towards Ising or XY). The peak value 
of C(T) is only weakly dependent on anisotropy, espe- 
cially for the Isinglike systems. A strong contrast to this 
behavior will be seen, as another quantum mechanical 
phenomenon, in the ferromagnetic XXZ chain. 

The linearity, at low temperatures, of the spin liq- 
uid (|i?| < 1) specific heat with respect to temperature 
is expected on the basis of spin-wave calculations for 
the antiferromagnetic XXZ model [45|, [4(|. This linear 
form of C(T) reflects the linear energy-momentum dis- 
persion of the low-lying excitations, the magnons. The 
low-tempcrature magnon dispersion relation is hu = ck n , 
where c is the spin-wave velocity and n = 1 for the an- 
tiferromagnetic XXZ model in d = 1 [40]. The internal 
energy, given by U = (1/2tt) / °° dkhw(k) / (e 15 ' 1 ^ -1), is 
dominated by the magnons at low temperatures, yielding 
U ~ T 2 and C ~ T for n = 1 in the dispersion relation. 
From this relation, our calculated spin-wave velocity c as 
a function of anisotropy R is given in Fig and compares 
well with the also shown exact result [47j . A simultaneous 
fit to the dispersion relation exponent n, expected to be 
1, yields 1.00 ± 0.02. However, for the Isinglike —R > 1, 
the unexpected linearity instead of an exponential form 
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FIG. 4: Calculated specific heats C of the antiferromagnetic 
XXZ chain, as a function of temperature for anisotropy coeffi- 
cients R = 0, -0.25, -0.50, -0.75, -1, -2, -4, -8, -oo span- 
ning the spin-liquid (upper panel) and Isinglike (lower panel) 
regions. 



0.50 



O 0.25 

+-> 

a 

| 0.00 
a 

CO 



0.25 



0.00 



0.0 



0.5 



1/ * Jxy 



0.0 0.5 0.0 0.5 



1.0 



1 

/X o 


i 

/ r X -0.5 

r \ 
■ / 

' / 

' m 

i / 
y 


1 

/^X 1 

ft t X 

' fx X. 
mi 

it 


i 

/\ 2 

// *\ 

■ I V X 

'/ 

1 / 



0.0 

V\J, 



0.5 



1.0 



FIG. 5: Comparison of our antiferromagnetic specific heat re- 
sults (thick lines) with the results of Refs.JH] (open circles), 
[f| (dotted), pj] (thin lines), [HJ (dash-dotted), and [USj] 
(dashed), for anisotropy coefficients R — 0, —0.5, —1, —2 span- 
ning the spin-liquid and Isinglike regions. 



caused by a gap in the excitation spectrum, points to the 
approximate nature of our renormalization-group calcu- 
lation. The correct exponential form is obtained in the 
large — R limit, where the renormalization-group calcu- 
lation becomes exact. 

Rojas et al. have obtained the high-temperature 
expansion of the free energy of the XXZ chain to order 
P , where (3 is the inverse temperature. The specific heat 
from this expansion is 



r = 2 + R 72 3R 3 

16 f " 32 xy 



6-8R 2 - R 4 



256 



-J, 



xy 



(17) 



This high-temperature specific heat result is also com- 
pared with our results, in Fig[6] and very good agreement 
is seen. In fact, when in the high-temperature region of 
< j3 < 0.1, we fit our numerical results for C(/3) to the 
fourth degree polynomial C — T,j =0 Ai(3' 1 , and we do find 

(1) the vanishing A < 10~ 5 and A\ < 1CT 7 for all R and 

(2) the comparison in Fig[5] between our results for A-i 
and A 3 and those of Eq. fTT)) from Ref.[25[, thus obtaining 
excellent agreement for all regions of the model. 



CD 
O 

'u 

CP 
CO 



J 



■H-J 



0.50 
0.25 
0.00 
0.25 
0.00 



1 2 1 2 1 2 3 







1 1 

y -0.75 


r 

J -i 


i i / 

J 2 


1 1 



1 2 1 2 1 2 3 

\Jz\ 

FIG. 6: Comparison of our antiferromagnetic specihc heat 
results (thick lines) with the high-temperature J —* behav- 
iors (thin lines) obtained from series expansion in Ref . [25l] . 
for anisotropy coefficients R — 0, — 0.50, — 0.75, — 1, — 2, — oo 
spanning the spin-liquid and Isinglike regions. 



V. FERROMAGNETIC XXZ CHAIN 

For the ferromagnetic (i.e., R > 0) systems in Fig|T] 
the (sfs|) expectation value becomes rapidly negative 
at lower temperatures for R < 1, even though for 
R > all couplings in the Hamiltonian are ferromag- 



netic. This is actually a real physical effect, not a nu- 
merical anomaly. In fact, we know the spin-spin corre- 
lations for the ground state of the one-dimensional XY 
model (the R = case of our Hamiltonian), and we can 
compare our low-temperature results with these exact 
values. The ground-state properties of the spin-i XY 
model are studied by making a Jordan- Wigner trans- 
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FIG. 7: Our calculated antiferromagnetic spin-wave veloc- 
ity c versus the anisotropy coefficient R. The dashed line, 
27rsin(7)/7 where 7 = cos _1 (— R), is the exact result [47| . 
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FIG. 8: Comparison of our results with the high-temperature 
expansion of Ref. 25| for all regions: antiferromagnetic (outer 
panels) and ferromagnetic (inner panels), spin-liquid (left 
panels) and Isinglike (right panels). Triangles and circles de- 
note our results, while solid and dashed lines denote the re- 
sults of Ref. 25] for A2 and A3, respectively. The error bars, 
due to the statistical fitting procedure of the coefficients A2 
and A3, have half-heights of 1.7 x 10 -4 and 2.6 x 10 -3 respec- 
tively. 



formation, yielding a theory of non-interacting spinlcss 
fcrmions. Analysis of this theory yields the exact zero- 
temperature nearest-neighbor spin-spin correlations [H 
shown in Table IIIII Our renormalization-group results 
in the zero-temperature limit, also shown in this table, 
compare quite well with the exact results, as with the 
other exact points in Fig. 1(b), although in the worst re- 
gion for our approximation. Finally, by continuity, it is 
reasonable that for a range of R positive but less than 
one, the z component correlation function is as we find, 
intriguingly but correctly negative at low temperatures. 
Thus, the interaction sfs* + sfs^ (irrespective of its sign, 
due to the symmetry mentioned at the end of Sec.IIB) 
induces an antiferromagnetic correlation in the s| com- 
ponent, competing with the s?s| interaction when the 
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2 4 6 0.0 0.4 0.8 

l/J z 

FIG. 9: Calculated nearest-neighbor spin-spin correlations 
(s^s^) (upper panels) and (sfsj) (lower panels) for the 
ferromagnetic XXZ chain, as a function of temperature, 
for anisotropy coefficients R = 0,0.25,0.50,0.75,1,2,4,8,00 
spanning the spin-liquid (left panels) and Isinglike (right pan- 
els) regions. 




l/J xy R 

FIG. 10: Left panel: Calculated nearest-neighbor spin- 
spin correlations (s|s|) for the ferromagnetic XXZ chain, 
as a function of temperature 1/J xy in the spin liquid, for 
anisotropy coefficients R = 0,0.25,0.50,0.75, 1. Right panel: 
The sign-reversal temperature To of the nearest-neighbor cor- 
relation (sfsj): our results (full curve) and the analytical re- 
sult from the quantum transfer matrix method (dashed) [23| ]. 



latter is ferromagnetic. 

For finite temperatures, our calculated nearest- 
neighbor spin-spin correlations are shown in Figs 03 1 1 01 
for different values of R. These results are compared with 
Green's function calculations [22j in FigQT] As expected 
from the discussion at the beginning of this section, in 
the spin- liquid region, the correlation (s?s|) is negative 
at low temperatures. Thus, a competition occurs in the 
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FIG. 11: Comparison of our ferromagnetic R = 1, | results 
with Green's function calculations [22] . 
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FIG. 12: Calculated specific heats C of the ferromagnetic 
XXZ chain, as a function of temperature for anisotropy coef- 
ficients R — 0, 0.25, 0.50, 0.75, 1, 2, 4, 8, 00 spanning the spin- 
liquid (upper panel) and Isinglike (lower panel) regions. 



correlation (sf s|) between the XY-induced antiferromag- 
netism and the ferromagnetism due to the direct cou- 
pling between the s z spin components. In fact, the rein- 
forcement of antiferromagnetic correlations of (sfs*) by 
increasing J xy (and also its converse) was seen in the 
antiferromagnetic XXZ chain discussed in the previous 
section. Thus, we see that whereas this cross-component 
effect is dominant at low temperatures in the ferromag- 
netic XXZ chain, it is seen at higher temperatures in the 
antiferromagnetic XXZ chain and, in-between, through- 
out the temperature range in the XY chain. 

In the ferromagnetic XXZ chain, as a consequence 
of the competition mentioned above, a sign reversal in 
(sfSj) occurs from negative to positive correlation, at 
temperatures J~j — Tn(R) .[4S| At this temperature, by 
cancelation of the competing effects, the nearest-neighbor 
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FIG. 13: Comparison of our ferromagnetic specific heat re- 
sults (thick lines) with the high-temperature J — > behaviors 
(thin lines) obtained from series expansion [25j], for anisotropy 
coefficients R = 0.25, 0.50, 0.75, 1, 2, 00 spanning the spin- 
liquid and Isinglike regions. 



Zero-temperature 
correlations of the 
spin-i XY chain 


Exact values 
from Ref. H 


Our 
RGT 
results 




0.15915 


0.17678 




-0.10132 


-0.12500 



TABLE III: Zero-temperature nearest-neighbor correlations 
of the spin- j XY chain. 



correlation (sfs|) is zero. Our calculated Xb(i?) curve is 
shown in Fig[T0l and has very good agreement with the 
exact result T = (\/3 sin.7/47) tan[7r(7r — 7V27] where 

7 = cos-^-i?) m. 

The calculated ferromagnetic specific heats are shown 
in Figll2l for various anisotropy coefficients and com- 
pared, in FigsH31 EH with finite-lattice expansion [||, 
quantum decimation I2l |. decoupled Green's functions 
[22| . transfer matrix |23l. |24j, high-temperature se- 
ries expansion (2f| results and, for the R = case, 
namely the XY model, with the exact result [5] C — 
(l/47rT) C (cosw/ cosh (^f^)) 2 dtu. In sharp contrast 
to the antiferromagnetic case in Sec. IV, the peak C(T) 
temperature is highest for the most anisotropic cases 
(XY or Ising) and decreases with anisotropy decreas- 
ing from either direction (towards Heisenberg). In the 
same contrast, the peak value of C(T) is dependent 
on anisotropy, decreasing, eventually to a flat curve, as 
anisotropy is decreased. This contrast between the fer- 
romagnetic and antiferromagnetic systems is a purely 
quantum phenomenon. Specifically, the marked contrast 
between the specific heats of the isotropic antiferromag- 
netic and ferromagnetic systems, seen in the full curves 
of Figsl4l and [T2l respectively, translates into the different 
critical t emp eratures of the respective three-dimensional 
systems. [UH^, H3, H^] Classical ferromagnetic and anti- 
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FIG. 14: Comparison of our ferromagnetic specific heat re- 
sults (thick lines) with the results of Refs. fal (d ash-double- 
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cles), and [23L I24I] (dashed), for anisotropy coefficients R — 
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ferromagnetic systems are, on the other hand, identically 
mapped onto each other. 

The low-temperature specifics heats are discussed in 
detail and compared to other results in Sec. VI. 
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FIG. 15: The calculated excitation spectrum gap A versus 
anisotropy. 
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VI. LOW-TEMPERATURE SPECIFIC HEATS 



FIG. 16: Calculated specific heat coefficient C/T as a function 
of anisotropy 7?, for T = 0.10, 0.05, 10~ 10 . 



Properties of the low-temperature specific heat of 
the ferromagnetic XXZ chain have been derived from 
the thermodynamic Bethe-ansatz equations [40]. For 
anisotropy coefficient \R\ < 1, the model is gapless 
fill [I2I ] and, except at R = 1, the specific heat is lin- 
ear in T — J~y in the zero-temperature limit, C/T = 
27/(3 sin 7) where again 7 = cos _1 (— R). Note that 
this result contradicts the spin-wave theory prediction of 
C ~ T 1 / 2 for the ferromagnetic chain (n = 2 for the fer- 
romagnetic magnon dispersion relation of the kind given 
above in Sec. IV). The spin- wave result is valid only for 
R = 1, the isotropic Heisenberg case. From the expres- 
sion given above, we see that C/T diverges as R — > 1~, 
and at exactly R = 1 it has been shown that C ~ T 1 / 2 

In the Isinglike region R > 1, the system exhibits a 



gap in its excitation spectrum and the specific heat be- 
haves as C ~ T _a exp(— A/T), with A being the excita- 
tion spectrum gap [TT|, [l^, ■ There exist two gaps for 
the energy, called the spinon ga p and the spin-wave gap, 
given by A spinm = \\/\- R~ 2 and A spinwave = 1-R' 1 . 
These are the minimal energies of elementary excita- 
tions [l(| S3| • A crossover between them occurs at R = | : 
below this value, the spinon gap is lower, while above 
this value the spin- wave gap is lower. We have double- 
fitted our calculated specific heats with respect to the 
gap A and the leading exponent a, for the entire range 



of anisotropy R between < R 



< 



1 (Figflgp. Our 



calculated gap A behaves linearly in for R^ 1 close 
to 1, and crosses over to 1/2 at R^ 1 = 0, as expected. 
We also obtain the exponent a = 1.99 ± 0.02 in the Ising 
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FIG. 17: Calculated specific heat coefficient C/T as a func- 
tion of temperature for anisotropy coefficient R — — 5 (thin 
grey), R = —2 (thick grey), -1 (dotted), -0.5 (dash-dotted), 
0.5 (dashed), and 2 (thin black). 



limit R- 1 < 0.2 and a = 1.52 ± 0.10 in the Heisenberg 
limit > 0.9. These exponent values are respectively 
expected to be 2 and 1.5 [9. Hoj. 

We now turn to the discussion of our specific heat re- 
sults for the entire ferromagnetic and antiferromagnetic 
ranges. Our calculated C/T curves are plotted as a func- 
tion of anisotropy and temperature in FigsfTB] and [T7I 
respectively. We discuss each region of the anisotropy R 
separately: 

(i) R > 1 : The specific heat coefficient C/T vanishes in 
the T — * limit and has the expected exponential form as 
discussed above in this section. The spin-wave to spinon 
excitation gap crossover is obtained. 

(ii) R w 1 : The double-peak structure of C/T in FigQJ] 
is centered at R = 1. As temperature goes to zero, the 
peaks narrow and diverge. 

(iii) — 1 < R < 1 : The specific heat coefficient is 
C/T = 27/(3 sin 7) in this region [Tl|,|4C]], and our calcu- 
lated specific heat is indeed linear at low temperatures. 
The C/T curves for R = -1,-0.5,0.5 in FigfT7] all ex- 
trapolate to nonzero limits at T = 0. The spin-wave 
dispersion relation exponent and velocity, for the antifer- 
romagnetic system, is correctly obtained for the isotropic 



case and for all anisotropies, as seen in Fig|7] FigfTSl 
directly compares C/T — 27/(3 sin 7) with our results: 
The curves have the same basic form, gradually rising 
from R = —1, with a sharp divergence as R nears 1. 
At R = 1+, we expect C/T = 0. Our T = 10~ 10 
curve diverges at R = 1 and indeed returns to zero at 
R = 1.0000001. 

(iv) R < — 1 : We expect a vanishing C/T, which we do 
find as seen in Fig|T5] and in the insets of Fig[T7] The 
exponential behavior of the specific heat is clearly seen 
in the Ising limit. 



VII. CONCLUSION 

A detailed global renormalization-group solution of 
the XXZ Heisenberg chain, for all temperatures and 
anisotropies, for both ferromagnetic and antiferromag- 
netic couplings, has been obtained. In the spin- liquid 
region, the linear low-temperature specific heat and, for 
the antiferromagnetic chain, the spin- wave dispersion re- 
lation exponent n and velocity c have been obtained. In 
the Isinglike region, the spin-wave to spinon crossover of 
the excitation spectrum gap of the ferromagnetic chain 
has been obtained from the exponential specific heat, 
as well as the correct leading algebraic behaviors in the 
Heisenberg and Ising limits. Purely quantum mechani- 
cal effects have been seen: We find that the xy corre- 
lations and the antiferromagnetic z correlations mutu- 
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FIG. 18: Calculated specific heat coefficient C/T as a func- 
tion of anisotropy coefficient R in the spin-liquid region, 
— 1 < R < 1, at constant temperature T = 10 -10 . Our 
renormalization-group result (grey curve) is compared to the 
zero-temperature Bethe-Ansatz result (black curve). Inset: 
our calculation (grey curve) at constant T = 10 -2 is again 
compared to the zero-temperature Bethe-Ansatz result (black 
curve) . 
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ally reinforce each other, for different ranges of temper- 
atures and anisotropics, in ferromagnetic, antifcrromag- 
netic, and XY systems. The behaviors, with respect to 
anisotropy, of the specific heat peak values and locations 
are opposite in the ferromagnetic and antiferromagnetic 
systems. The sharp contrast found in the specific heats 
of the isotropic ferromagnetic and antiferromagnetic sys- 
tems is a harbinger of the different critical temperatures 
in the respective three-dimensional systems. When com- 
pared with existing calculations in the various regions of 
the global model, good quantitative agreement is seen. 
Even at zero temperature, where our approximation is at 
its worst, good quantitative agreement is seen with exact 
data points for the correlation functions (FigQJb)), which 
we extend to all values of the anisotropy. Finally, the rel- 



ative ease with which the Suzuki- Takano decimation pro- 
cedure is globally and quantitatively implemented should 
be noted. 
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